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Abstract 

This paper is devoted to study the dynamics of gravitational col- 
lapse in the Misner and Sharp formalism. We take non- viscous heat 
conducting charged anisotropic fluid as a collapsing matter with cylin- 
drical symmetry. The dynamical equations are derived and coupled 
with the transport equation for heat flux obtained from the Mitller- 
Israel-Stewart causal thermodynamic theory. We discuss the role of 
anisotropy, electric charge and radial heat flux over the dynamics of 
the collapse with the help of coupled equation. 
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1 Introduction 



One of the most important problems in the gravitation theory and relativis- 
tic astrophysics is to understand the end state of a continual gravitational 
collapse. A massive star undergoes to gravitational collapse at the end of 
its life cycle. This happens when all the internal nuclear forces fail to sup- 
ply the sufficiently high pressure to counter-balance gravity. The compact 
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objects such as white dwarfs, neutron stars and black hole are the results 
of possible stages of the collapsing astronomical objects. In white dwarfs 
and neutron stars, gravity is neutralized by electron and neutron degeneracy 
pressure respectively and black hole is a complete collapsed object. 

Oppenheimer and Snyder (1939) are the pioneers who studied gravita- 
tional collapse of an adiabatically flowing dust. This was idealized problem 
because dust is unrealistic matter and one cannot ignore the effects of pres- 
sure on the formation of spacetime singularity. A more analytic analysis was 
made by Misner and Sharp (1964) with perfect fluid in the inner region of 
a star. They formulated the dynamical equations governing adiabatic rela- 
tivists collapse. In both cases, vacuum was taken in the exterior region of a 
star. 

The concept of non-vacuum exterior of a star was introduced by Vaidya 
(1951) for the radiating fluid in the interior region of the star. Goswami 
(2007) formulated a more realistic collapsing model by taking the radiat- 
ing dust matter in the interior of a star. He remarked that bounce in the 
collapse is due to the dissipation. Debnath et al. (2005) explored gravita- 
tional collapse of the non-adiabatic fluid by assuming quasi-spherical Szekeres 
spacetime in the interior and plane symmetric Vaidya solution in the exte- 
rior region. By using the local conservation of momentum, they studied the 
thermodynamical behavior of the collapsing matter. 

An extensive literature survey (Herrera et al. 2004 and Mitra 2006)pre- 
dicts that gravitational collapse is highly dissipative process. This indicates 
that the effects of the dissipation must be included in the study of collapse 
for its better understanding. Herrera and Santos (2004) explored dynamical 
description of gravitational collapse by using Misner and Sharp's formula- 
tion. Matter under consideration was distributed with spherically symmetric 
and energy loss in the form of heat flow and radiation. Chan (2001) studied 
the realistic model of radiating star which undergoes dissipation in the form 
of radial heat flow and shear viscosity. Herrera et al. (2009) also formulated 
the dynamical equations by including dissipation in the form of heat flow, 
radiation, shear and bulk viscosity and then coupled with causal transport 
equations with spherical symmetry. Herrera (2006) discussed the inertia of 
heat and its role in the dynamics of dissipative collapse with outgoing radial 
heat flux by using spherical symmetry. 

Most of the work available in spherical symmetry is due its simplest sym- 
metry. To generalize the geometry of the star, people worked on gravitational 
collapse using the non-spherical symmetry. The existence of cylindrical and 
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plane gravitational waves provides strong motivation in this regard. Herrera 
et al. (2005) formulated the set of equations with regularity and match- 
ing conditions for the static cylindrically symmetric distribution of matter. 
Sharif and Ahmad (2007) studied cylindrically symmetric gravitational col- 
lapse of two perfect fluids using the high speed approximation scheme. They 
investigated the emission of gravitational radiations from cylindrically sym- 
metric gravitational collapse. Nolan (2002) investigated naked singularities 
in the cylinderical gravitational collapse of counter rotating dust shell. Di 
Prisco et al. (2009) discussed the shear free cylindrical gravitational collapse 
by using junction conditions. Nakao et al. (2009) studied gravitational col- 
lapse of a hollow cylinder composed of dust. Recently Sharif and Rehmat 
(2010) discussed the dynamics of viscous dissipative plane symmetric gravi- 
tational collapse. 

The behavior of electromagnetic field in gravitational field has been the 
subject of interest for many people. Thorne (1965) developed the concept 
of cylindrical energy and investigated that a strong magnetic field along the 
symmetry axis may halt the cylindrical collapse of a finite cylinder before 
it reached to singularity. In recent papers (Sharif and Abbas 2009, 2010a, 
2010b), we have studied the effects of the electromagnetic field on the gravi- 
tational collapse by taking the homogenous, non-homogeneous and spherical 
model. Di Prisco et al. (2007) derived the dynamical equations for the spher- 
ically symmetric collapse by including electromagnetic field. This work has 
been extended by Sharif and Siddiqa (2011) for the charged plane symmetric 
gravitational collapse. Also, Sharif and Fatima (2011) discussed dynamics of 
adiabatic charged viscous cylindrical gravitational collapse. 

This paper is aimed to study the dynamics of non-adiabatic charged cylin- 
drically symmetric gravitational collapse to see the effects of charge and heat 
flux on the process of collapse. The plan of the paper is the following. In the 
next section, we describe the gravitational source and the Einstein-Maxwell 
field equations. Section 3 is devoted to matching conditions. We formulate 
the dynamical equations in section 4 and the derivation of the transport 
equation and their coupling with the dynamical equations are presented in 
section 5. The last section contains the conclusion of the paper. 
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2 Interior Matter Distribution and the Field 
Equations 

We take non-static cylindrically symmetric as an interior metric in the co- 
moving coordinates in the form 

ds l = _A 2 dt 2 + B 2 dr 2 + C 2 (d9 2 + dz 2 ), (1) 
— oo < t < oo, ^ r < oo, 0^6^ 2%, — oo < z < +oo 

where A, B and C are functions of t and r. Matter under consideration 
is anisotropic fluid which undergoes dissipation in the form of heat flux. 
The energy-momentum tensor for such a fluid dissipating only at diffusion 
approximation, i.e., e = is defined as (Herrera 2006)) 

T af} = (fi + P ± )V a Vp + P ± g a p + (P r - P ± )x a X/3 + V a qp + Vpq a , (2) 

where fx, P r , P±, q a , V a and \ a are the energy density, the radial pressure, 
the tangential pressure, heat flux, the four-velocity of the fluid and the unit 
four- vector along the radial direction respectively. For the metric (1), the 
four-vector velocity, heat flux and unit four-vector along the radial direction 
are given by 

V a = A-X, X a = B- l 5 a ^ q a = B- l q5^ 

which satisfy 

v a v a = -i, x a x a = i, x a v a = o, q a V a = 0. 

We can write the electromagnetic energy-momentum tensor in the form 

= ^ (Wr - J) ■ (3) 

The Maxwell equations are given by 

F a/ 3 = 4>/3 >a ~ 4>a,p, (4) 

F a % = 4vrJ a , (5) 

where F a/3 is the Maxwell field tensor, (fi a is the four potential and J a is 
the four current. Since the charge is at rest with respect to the co-moving 



4 



coordinate system, thus the magnetic field is zero. Consequently, the four 
potential and the four current will become 

= K> ja = ° va > ( 6 ) 

where (ft = r ) is an arbitrary function and a = a(t, r) is the charge 
density. 

For the interior spacetime, using Eq.flB]), the Maxwell field equations take 
the following form 

/ a' w m\ 

= AnaAB 2 , (7) 

t = 0, (8) 

where dot and prime represent the partial derivatives with respect to t and 
r respectively. Integration of Eq.([7j) implies that 
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where s (r) = 2ir J^aBC 2 dr is the total charge distributed per unit length of 
the cylinder and is the consequence of law of conservation of charge, J!i = 0. 
Obviously Eq.flH]) is identically satisfied by Eq.(j9]). 

The Einstein field equations, G a $ = 87r(T a/3 + T^™^), for the metric (1) 
can be written as 

84T 00 + T ( em) ) = 8ttM 2 + ^ 

C ( B c\ fA\ 2 f C" B'C ,C' 

= c[ 2 b + c) + {b) (- 2 c +2 bc- ( c) 

(10) 

/m m ( Pm \, 2 / C BC CA' \ 

8.(T 01 + Tr ) = Sn q = — - _ - _ j , (11) 

B\ 2 ( C ( cY AC\ fC'V A'C 
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87r(T 22 + E 22 ) = 8ttP ± C 2 + ^ 



C\ 2 I B C A I B C\ BC 



A) \B + C A\B + C I ~ BC 

C\ 2 (A" C" _A' (B 1 C'\ B'C'\ 
Bj \A + ~C'H\B~~C)~^C)- [1?>) 

The C-energy for the cylindrically symmetric spacetime is defined by (Thorne 
1965) 

E=l(l-r 2 V a rV a r), (14) 

o 

The circumference radius p, specific length / and areal radius f can be defined 

as 

P 2 = £(l)a£fl); I 2 = ^(2)a^(2), r — pi, 

where £(i) = ^ 2 ) = J| are Killing vectors and E represents the gravita- 
tional energy per unit specific length of the cylinder. 

The specific energy of the cylinder (Poisson 2004) analogous to Misner 
and Sharp energy for the spherical symmetry in the interior region with the 
contribution of electromagnetic field can be written as follows 

„ ( etc 2 C a \ s 2 

E= 8 + l[lp-&) + 2C- (15) 

We would like to mention here that this energy is also analogous to Taub's 
mass function in the plane symmetric spacetime (Sharif and Rehmat 2010 ). 



3 Junction Conditions 

In this section, we assume that the 3D timelike boundary surface £ splits 
the two AD cylindrically symmetric spacetimes V + and V~. The metric 
which describes the internal region V~ is given by Eq.([T]) while for the rep- 
resentation of exterior region V + , a metric in the retarded time coordinate 
is considered. If M{y) and Q{v) are mass and charge, respectively, in re- 
tarded time then the corresponding cylindrically symmetric spacetime given 
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by (Chao-Guang 1995) will take the form 

ds l = _ ( ~2M(t/) + Q^± )d]y 2 _ 2dRdu + R 2 {d0 2 + dz 2^ (;|G) 

ix R 



where v is the retarded time coordinate. We can write the induced metric 
for the hypersurface £ in the following form 

(ds% = -dr 2 + y 2 (r)(d6 2 + dz 2 ), (17) 

where = (r, <fi, z) (i — 0, 2, 3) represent the intrinsic coordinates of E. 

The Darmois junction conditions (Darmois, 1927) can be stated as fol- 
lows: 

• The continuity of the first fundamental form over the hypersurface £ 
i.e., 

(ds 2 )x = (ds 2 _)x = (18) 

• The continuity of the second fundamental form over the hypersurface 
£ 

[tfy] = K± -K- J= 0. (19) 
Here, Kf- is the extrinsic curvature given by 

a2 a d\ v 

x 5 = -^ ( W + r -## ) ' ^^- = ' 1 ' 2 ' 3 )- ( 20 ) 

where n^r are the components of outward unit normal to the hypersurface in 
the coordinates x ±M - 

We can write the equations of hypersurface as follows 

h-(t,r) = r-r s = 0, (21) 
h + (u,R) = R-Rx(v)=0, (22) 

where rs is a constant. Using above equations, we have the interior and 
exterior spacetimes on £ as follows 

(cfe 2 _) s = -A 2 (t,r^)dt 2 + C 2 {t,r^)(d6 2 + dR^)2). (23) 
+ Rl(d9 2 + dz 2 ). (24) 
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The continuity of the first fundamental form gives 

R s (u) = C(t,r s ), (25) 

dt 1 . , 

Tr = A' (26) 
dv u -2M(y) Q», 2dR^ 

Tt ~ K— R^ + _ RT ) + ^r J ■ (27) 

Now we consider the second fundamental form over E. For this purpose, we 
need the outward unit normals to £ using Eqs. (l21j) and 



n~ = 5(0,1,0,0), (28) 
-2M{v) Q 2 (v) 2dR^\^ , dR^ 



< = + + (--sr- 1 - - )- (29) 



The non-zero components of the extrinsic curvature Kf- are 
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The continuity of the extrinsic curvature components yields 

Using Eqs.fl25J-(J27]), (HU) and §TQ in Eqs.flS]) an d it follows that 

E'-M = -os = Q, (36) 
8 

I^Pr-g. (37) 
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The first equation indicates that the difference between two masses is equal 
to | as shown in the adiabatic case (Sharif and Fatima 2011). This is due to 
the least unsatisfactory definition of Thorne C-energy (Thorne 1965). The 
second equation describes a relation between heat flux, radial pressure and 
charge over the hypersurface E. It is obvious from this equation that for 
uncharged radiating fluid, radial pressure and heat flux are equal over the 
boundary of the collapsing cylinder. 



4 The Dynamical Equations 

Here we derive the dynamical equations for non-adiabatic charged anisotropic 
fluid. The energy-momentum conservation, (T a/3 + T^ em ) a ^) ;/ 3 = 0, implies 
that 
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(39) 



BC y ~' tiBC a 

Following Misner and Sharp formalism (Misner and Sharp 1964), we discuss 
the dynamics of the collapsing system. We introduce proper time derivative 
as well as the proper radial derivative constructed from the circumference 
radius of a cylinder inside E as (Sharif and Fatima 2011) 

l__d l_d_ 

Adi' R ~~R'!h' 
The fluid velocity in the case of collapse can be defined as 



D 



T 



R = C. 



(40) 



U = D T (R) = D T (C) < 
Consequently, we can write Eq. (jT5!) as 

s 2 2 / , I 
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Using Eqs. flTT]) . ( 1T2|) . (lT5"j) and (j4"0~|) . the time rate of change of the C-energy 
turns out to be 



D T E' 



-AttR 2 



P 



1 



32nR 2 



U + qE 



3s 2 U 
2R 2 ' 



(43) 



This equation represents the variation of total energy inside the collapsing 
cylinder. Since U < 0, the first term on the right hand side of this equation 
will increase the energy of the system provided that the factor within the 
round brackets is positive. The second term in the square brackets due to 
negative sign describes the outflow of energy in the form of radiation during 
the collapse. For the collapsing cylinder containing the same species of the 
charges, the third term will decrease the energy of the system as plays 
the role of Coulomb repulsive force and U < 0. 

Similarly, using Eqs.fllO]), ([II]), (JT5J) and fllO]), we obtain 



D R E' 



U. I s „ 3s 
D R s + 



AnR 2 (fi + q^) + 
E 



8 + R' 



2R 2 ' 



(44) 



This equation gives the variation of energy between the adjacent cylindrical 
surfaces inside the matter distribution. The first term is the energy density 
of the fluid element along with heat flux contribution. Since U < 0, the heat 
flux factor decreases the energy of the system during the collapse of cylinder. 
The term | comes from the definition of C-energy and the remaining terms 
are due to the electromagnetic field. 

Using Eqs.( TT2j) . (JT5]) . (14"T|) and we can obtain the acceleration D T U 
of a collapsing matter inside E 



_ EA' 5s 2 
AnRP r + — — + 



AB 2R 3 ' 



(45) 



Inserting the value of 4^ from Eq. (j43|) into Eq. (l39|) . it follows that 
(/j + P r )D T U= - (m + P, 
- E 



R 2K 8 J 2R 3 



AqU 
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(46) 
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where G = — %) ^ for simplicity. Now the complete dynamics of 
the system is described by Eg. (|46p . The system will evolve radially outward 
or inward according to DjU < or DjU > 0. Thus the terms in Eq. (j46j) 
contributing negatively, favors the collapse while the other contribution pre- 
vents the collapse. If both of these cancel each other then there will be a 
hydrostatic equilibrium. Since the left hand side of Eq.( H6|) represents force, 
so the factor /i + P r refers to an inertial mass density independent of charge 
and heat flux contributions. The first and third terms on the right hand side 
represents the gravitational force. The second term represents the heat flux 
contribution which seems to leave the system (due to negative sign) through 
the outward radially directed streamlines. Being in the same direction of 
pressure, it supports the pressure and would prevent the collapse. 

The term (fi + P r ) [-^(E' — |) + AnP r R — |^] represents the gravitational 
force. The factor within the first square brackets shows the effects of specific 
length and the electric charge on the active gravitational mass term (fi + P r ). 
The third term has three main contributions, i.e., the first is the pressure 
gradient which is negative, the second is the local anisotropy of the fluid 
which will be negative for P r < P^ and the third is the electromagnetic field 
term. For an isotropic pressure, the second contribution will be vanished. 
Further, following Di Prisco et al. (2007), it can be found that the third term 
contributes negatively for -| > Drs. Thus the third square brackets, under 
these conditions with negative sign, contributes positively by reducing an 
attractive behavior of force appearing on left hand side of the equation. Since 
the attractive force is decreased, so the third term prevents the gravitational 
collapse of the cylinder. 

5 The Transport Equation 

The transport equation predicts the processes of mass, heat and momentum 
transfer during the dynamics of a realistic matter. The transport equation 
for heat flux derived from the Mitller-Israel-Stewart causal thermodynamic 
theory is given by (Herrera 2006) 

1 rV b 

rh ab V c q b , c + q a = -Kh a \T b + a b T) - - K T 2 (—). <b q a , (47) 

where h ab = g ab + V a V b is the projection tensor, k denotes thermal con- 
ductivity, T is temperature, r stands for relaxation time which is the time 
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taken by a perturbed system to return into an equilibrium state and a^T is 
the Tolman inertial term. Due to symmetry of the spacetime, the transport 
equation reduces to the following form 

ri = -r^y - i, 9 (| + 2§) - ±(TAY - q A. (48) 



Using Eqs.f HOj) and (jUj) in this equation, it follows that 

^ nT 2 q „ . r , r 3U _ 1, kF „ m 

x ^-^ + 1+4^-^. (49) 

In order to understand the effects of heat flux or dissipation on collapsing 
process, we couple Eq.(f49l) with dynamical Eq.()46l). Thus the replacement 
of Eq.pU in Eq.gED, yields 

(fi + P r )(l-a)D T U = (l- a )F grav + F hyd + ^D R T 

+ ^*<^>l " + G " ?)■ < 50 > 

where F sra „, and a are given by the following equations 

= -(^ + P r )[i?'-^+47rP r i? 3 -^]-l, (51) 

= -E 2 [D R (P r ) + ^(Pr ~ P±) ~ ^^], (52) 
a = — (fi + Pr)- 1 . (53) 

T 

For the physical interpretation of Eq.f l50p . it can be observed that the left 
hand side of this equation being the product of inertial mass density (/i + 
P r )(l — a) and acceleration, Dj^U represents the Newtonian force. Thus we 
can write F = (/i + P r )(l — q)DtU ■ It is clear that when a — > 1, then F — > 0, 
which means that there is no inertial force and matter would experience the 
gravitational attraction which causes the collapse. For < a < 1, the 
inertial mass density goes on decreasing while 1 < a indicates the increase 
of inertial mass density. Of course, by the equivalence principle, there would 
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occur decrease and increase in the gravitational mass. In this way, one can 
explicitly distinguish the expanding and collapsing mechanism during the 
dynamics of dissipative system. 

Also, Eq.( l50p implies that gravitational force is affected by the same fac- 
tor but hydrodynamical force is independent of this factor. Further, as long 
as (n + -P r )(l — oc)DtU < 0, the total Newtonian force of the system remains 
directed downward which is the indication for the gravitational collapse and 
the converse is true for the expansion. If there is a continuous change in a 
from a value greater than one to less than one and vice versa, then there 
is a transition phase in the system and bouncing would occur. This phe- 
nomenon causes the loss of energy from the system and hence the collapsing 
cylinder with non-adiabatic source leads to the emission of the gravitational 
radiations. On the basis of this fact, the exterior of the collapsing cylinder 
is radiation zone which is completely described by a spacetime in radiation 
coordinates like in the present case. 

6 Outlook 

This paper deals with the effects of the charge and heat conduction on the dy- 
namics of cylindrical anisotropic fluid collapse. We have extended the recent 
work of Sharif and Fatima (2011) to non-adiabatic case for the transporta- 
tion process of heat flux during the dynamics of realistic matter collapse. For 
this purpose, the non- viscous heat conducting anisotropic fluid with cylindri- 
cal symmetry has been taken as the source of gravitation in the presence of 
electromagnetic field. Using the Misner and Sharp formalism, the dynamical 
equations are derived. We have found that during the collapse of non-viscous 
heat conducting charged anisotropic fluid, the radial heat flux and electric 
charge causes to reduce the energy of the system. 

For q > 0, £" = ^->0, Eq. (|44p yields that the second term in the first 
square brackets is negative which indicates that heat is emitting from the 
system. Thus the external region of such a collapsing system, being non- 
vacuum, is defined by a cylindrical geometry in the radiation (single null) 
coordinate. This prediction is analogous to the Vaidya (1951) for spherically 
symmetric case. The time evolution of the system, given by dynamical equa- 
tion ( )4"6]) , indicates that the inertial mass density is independent of heat flux 
and electric charge. The left hand side of this equation corresponds to the 
Newtonian force of the system which is decreased by the heat flux. 
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Since the collapse of a star is an irreversible process. The transport pro- 
cess of such non-equilibrium objects and connection between their dynamics 
and thermodynamics are important for the better understanding of the prob- 
lem. Thus using the Miiller-Israel-Stewart causal thermodynamic theory, the 
transport equation for the dissipative fluid has been formulated and coupled 
to the dynamical equation. The coupled equation helps to determine the in- 
fluence of the heat flux over the dynamics of collapsing cylinder. It has been 
found that in the coupled dynamical Eq. (l50|) . the inertial as well as gravita- 
tional masses are influenced by the factor (1 — a). The role of a during the 
dynamics of system can be explained as follows: For a tends to one, we get 
zero mass density. For < a < 1, the inertial and the gravitational mass 
density are decreased while for a > 1, the gravitational force term becomes 
negative. This is the case for the reversal of collapse. The conditions on a 
have been investigated for the bouncing behavior of the system. 
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